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Abstract: Motivated by WMAP's confirmation of an anomalously low value of the 
quadrupole moment of the CMB temperature fluctuations, we investigate the effects on 
the CMB of cutting off the primordial power spectrum P(k) at low wave numbers. This 
could arise, for example, from a break in the inflaton potential, a prior period of matter 
or radiation domination, or an oscillating scalar field which couples to the inflaton. We 
reanalyze the full WMAP parameter space supplemented by a low-/c cutoff for P(k). The 
temperature correlations by themselves are better fit by a cutoff spectrum, but including 
the TE temperature-polarization spectrum reduces this preference to a lAa effect. Inclu- 
sion of large scale structure data does not change the conclusion. If taken seriously, the 
low-/c cutoff is correlated with optical depth so that reionization occurs even earlier than 
indicated by the WMAP analysis. 
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1. Introduction 

One of the intriguing results of the Wilkinson Microwave Anisotropy Probe (WMAP) is 
a smaller than expected correlation of temperature fluctuations on large angular scales, 
corresponding to a low value for the quadrupole moment. Figure 1 shows the WMAP 
results for C(9) and the low-multipole Q's. Their "Basic results" paper jjj states that the 
"probability of so little C{9 > 60°) anisotropy power is ~ 2 x 10~ 3 , given the best-fit [run- 
ning spectral index] ACDM model." (See also refs. [^, ||] for related discussions.) Given the 
large uncertainties in this region due to cosmic variance, one might never know whether this 
constitutes a truly significant deviation from standard cosmological expectations. However 
it has recently been suggested that measuring polarization from galaxy clusters may in the 
next decade give a complementary determination of the quadrupole moment 0] and its 
evolution with redshift. 

Loss of power at large angles could be explained by lowering the primordial power 
spectrum of inflaton fluctuations (P(k)) at small wave numbers. Prior to the WMAP ob- 
servation, there have been several suggestions for generating such an effect. Sharp features 
in the inflaton potential ||-@ can suppress power at low k. Recently |§ studied two other 
possibilities: an oscillating scalar field % coupling to the inflaton can temporariliy suppress 
inflaton fluctuations, leading to a cutoff at low k; a prior period of matter or radiation 
domination before the beginning of inflation has a similar effect ||. Spatial compactness 
|1C] or curvature [|ll]] can also lead to a smaller CMB quadrupole moment. These scenarios 



require that the total amount of inflation is close to the minimum amount required to solve 
the horizon and flatness problems, or alternatively that the break in the inflaton potential 
is reached just when the relevant scales are crossing the horizon. 



Recent papers [12, 1 3 1 considered how cutting off P(k) at low k can improve the fit to 



the data. Ref. [13] found that the data favor having such a modified primordial spectrum at 



the 2a level, whereas [12] obtain weaker evidence for a cutoff. The conclusion of ref. [|13j was 
based upon using only the temperature anisotropy data, and ignoring polarization. The 
TE cross power spectrum represents one third of the total data set, and has an important 
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Figure 1: Left: Angular correlation function C(8) for CMB temperature anisotropy, from fig. 13 
of Q. Our plot of WMAP TT multipole data (from http://lambda.gsfc.nasa.gov/product/map/ 
map_tt_powspec.cfm) with the best fit ACDM model (with constant scalar index n s ) from table 7 
of Jl5| . Cosmic variance is a factor of 1 ± 1/^/1 + 1/2. 



effect on the determination of parameters that would lower the TT quadrupole moment. 
If one simply cuts off P(k) at low k while keeping all other cosmological parameters fixed, 
the improvement of the fit to the TT spectrum is accompanied by a deterioration in the 
fit to the TE spectrum, such that the overall improvement of the fit is nil. To counteract 
this, one needs to allow other parameters like the optical depth r and the fraction of dark 
energy £l\ to vary. In the present study, we do a complete likelihood analysis in which all 
the relevant cosmological parameters are varied in order to find those values which have 
the real maximum likelihood. We corroborate and extend the results of [O], who used a 
similar approach. 

In the second section we describe the theoretical models which can lead to a suppression 
of low- A; power in the primordial fluctuations. Section three presents our likelihood analysis, 
in which we find that the preference for a cutoff in P(k) is of marginal statistical significance. 
Our findings and conclusions are summarized in the final section. 

2. Theoretical Models 

There are many ways to alter the spectrum of inflaton fluctuations relative to the flat 
Harrison-Zeldovich form. Here we will concentrate only on those that produce a sharp 
reduction in large wavelength power, needed to suppress the low multipoles of the CMB. 
One way is to engineer the inflaton potential V{4>), which in principle can yield any desired 
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form for P(k) |g, 14]. A model which predicts a step-like feature in P{k) was proposed 
by Starobinsky |J], which assumes that there is a sudden change in the slope of V (</>)• 
If by chance the scales presently corresponding to large-angle CMB anisotropies exit the 
Hubble radius at the moment when cp crosses the kink in V (</>), then the power of smaller k 
fluctuations which subsequently cross the horizon can be suppressed. The resulting power 
spectrum (shown in Jiffl) looks similar to the one we shall discuss below in fig. 2b. 

More recently (but before WMAP's data release) ref. M investigated possible effects 
of a limited duration of inflation on the power spectrum of the inflaton fluctuations. Two 
of these predicted that P(k) should be strongly suppressed at low wave numbers. The first 
was in the context of the hybrid inflation model, with Lagrangian 



C 



(2.1) 



with V(4>, x ) = \m 2 4> 2 + \ X( X 2 ~ v 2 ) 2 + \ g X V + ^ A </> 



It was noticed that if the field which is responsible for triggering the end of inflation, is 
oscillating prior to horizon crossing of the relevant inflaton fluctuations, these oscillations 
can strongly affect P(k) for k values which are below some characteristic scale k max , defined 
by 

U = Me- Ht0 (2.2) 

Here M is the \ field mass, and to is the amount of time prior to horizon crossing of the 
mode k during which x was oscillating. In fact there are two important scales, this one, 
which comes about because the x oscillations cause production of inflaton fluctuations, and 
a lower one, below which the fluctuations are suppressed. The suppression occurs because 
the average value of gx 2 contributes to the squared mass of the inflaton, dm^ = (gx 2 ) = 
^9Xo e ~ 3Ht ■ When Sm 2 ^ ^ H 2 , the inflaton rolls quickly, which suppresses fluctuations on 
scales k which have not yet crossed the horizon. However, the amplitude of x redshifts 
exponentially, and at a certain time Srrv^ will fall below H 2 . Fluctuations which have not 
yet frozen out by this time will have a chance to grow toward their normal amplitude, 
5(j)(k) = H/2tt; thus larger wave numbers will be unchanged by this effect. The fractional 
amount by which P(k) is reduced at low k is gx1\/H 2 f° r small values of this parameter, 
where xo is the initial amplitude of the x oscillations, and H is the Hubble parameter during 
inflation. For large values of gxo/H 2 , the suppression is essentially complete, P(k) = 0. 
This is illustrated for some typical parameter values in figure 2a. 

The second situation studied in || is that in which the inflaton fluctuations with wave 
number k cross the horizon shortly after the beginning of inflation, where the inflationary 
epoch is preceded by matter or radiation domination Q . In this case, power is suppressed 
on scales below k ~ H, the only relevant scale in the problem. (This refers to the value 
of k at the time when inflation starts. The physical wave number gets reduced by the 
subsequent inflation.) The deformation of the power spectrum has a unique form, which 
is shown in figure 2b for the case of prior matter domination. The suppression of P(k) in 
this case is a factor of 0.063. 

Both of the above effects rely upon having a limited period of inflation: inflation 
must not have started much earlier than the time when the large-angle CMB fluctuations 
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Figure 2: (a) Effect of x oscillations on P(k) for M = 50H, gxl/e 3Ht ° = 0.003H 2 , to = 4/F. 
P(k) is normalized such that P(k) = 1 for a flat Harrison-Zeldovich spectrum, (b) Effect of prior 
period of matter domination on P(k). The normalization of k is in fact arbitrary — see discussion 
at beginning of section 3. 



first crossed the inflationary horizon. For the oscillating scalar scenario, this is because 
its oscillations are quickly Hubble damped. After at most 30 e-foldings of inflation, the 
amplitude of x wm be too small to have any further effect on the CMB. In the second 
scenario, there must be a coincidence of scales such that the large angle fluctuations cross 
the horizon just at the beginning of inflation; in this case inflation can last no longer than 
the minimum duration needed for solving the horizon and flatness problems. 



3. Effect on CMB temperature fluctuations 

Let us consider the effect of the spectrum shown in fig. 2b on the CMB temperature 
anisotropy. The physical scale of wave numbers at which P(k) changes is determined by 
the total amount of inflation following the horizon crossing of the corresponding mode, 
denoted by k = H in fig. 2b. A longer period of inflation stretches the wavelength of this 
mode to larger physical values. We refer to this model-dependent cutoff scale as fc*. Figure 
3 shows the effect of varying k* in a modified version of the CMBFAST code |l6|] which 
incorporates the distorted spectrum. There it is seen that the relevant scales for suppressing 
the lowest multipoles are on the order of fc* ~ a few xlCT 4 Mpc" 1 . As noted in fTJ, @> 
even a sharp cutoff in P(k) does not generate a sharp rise in the temperature multipoles 
since the latter are a convolution of the former of the form C\ ~ f dkk~ l jf{k[r]Q—r]^ ec \)P{k). 

If we ignore the TE polarization data, it is possible to obtain a better fit to the 
measured WMAP anisotropy in the manner shown in fig. 3. However, if one changes P(k) 
like this without altering any other parameters, it exacerbates the fit to the TE spectrum. 



Figure 4a shows the experimental data at low I [ IT | ; fig. 4b shows that larger values of 
the cutoff wave number fc* suppress the low TE multipoles in conflict with experimental 
observation of stronger power at low I. The strong TE signal at low I is the basis for 
WMAP's inference of a large optical depth r and early reionization, which shows that r 
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Figure 3: Effect of prior period of matter domination on temperature anisotropy. The cutoff wave 
number fc* is in units of Mpc -1 . 



is an important parameter to vary in order to try to repair the damage to the fit to TE 
done by modifying P(k). However changing r can also hurt the agreement between data 
and theory for the higher multipoles unless the scalar spectral index n s is also adjusted to 
compensate the change in r. 
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Figure 4: (a) WMAP TE cross polarization spectrum (fig. 8 of 0). (b) Suppression of TE 
spectrum by P(k) with low-fc cutoff fc*. 

This kind of reasoning suggests that no safe conclusions can be drawn without doing a 
complete analysis in which at least a minimal set of cosmological parameters (like the six- 
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dimensional parameter space considered in []l], |l5|) are allowed to vary. We have therefore 
generated grids 1 of theoretical power spectra with varying values for cjj = f^/i 2 , uj m = 
r2 m /i 2 , h, r, A s , n s and for the cutoff scale fe*, using several forms for the cutoff spectrum 
P(k). We fixed to zero some parameters which are not strictly necessary for fitting the 
WMAP data, like the spatial curvature or the amount of primordial gravitational waves. 
In the following, the minimum x 2 values have been found using a Powell minimization 
algorithm. For models which do not match one of the grid points, our code first computes 
the power spectra Cj" 1 ', Cf E , P{k) by cubic interpolation between the grid's neighboring 
points, and then calculates the corresponding x 2 value. We checked that for our grid 
spacing, the interpolation is always accurate to one per cent. To be sure that we understand 
the results of ref. (T^j, we started with their ansatz 

P(k) = A s (l - e - ( fc/fc * )a )fc n '- 1 (3.1) 

where a = 3.35 is chosen to match the shape at low k of a model similar to that shown 
in fig. 2b (they considered a period of prior kination rather than matter domination). We 
performed comparisons with three different data sets: WMAP TT power spectrum alone, 
WMAP TT and TE combined spectra, and WMAP TT+TE supplemented by large scale 
structure data from the 2 degree field (2dF) galaxy redshift survey |l8|]. For WMAP, we 
computed the likelihood of a given model using the software provided by the collaboration 



and described in [IS]. For the 2dF power spectrum, we used the window functions and 
correlation matrix available at http://msowww.anu.edu.au/2dFGRS/. 

First, we verify that there is some preference for a nonvanishing cutoff in k space when 
only the temperature data is included. Figure 5 shows the likelihood (normalized to 1 
at the best fit value) for k*, having marginalized over the other parameters. The model 
without a cutoff is excluded only at the 90% confidence level (A% 2 ~ 2.6), which is not 
very significant. We note that artificially changing the first three theoretical multipoles, 
so as to exactly match the first three data points, improves x 2 by 7.0, corresponding to 
a Bayesian confidence level of more than 99%. Had we been able to achieve such a large 
Ax 2 , we would have been more convinced that there is some evidence for a cutoff power 
spectrum in the temperature data. 

We obtain a smaller preferred value of k* = 0.0003 Mpc" 1 than did ref. llDf ], who found 
k* = 0.0005 Mpc -1 . We have verified explicitly that this discrepancy is due to their having 
done the analysis on a smaller grid, in a restricted subspace of the parameters, which does 
not include the true minimum x 2 ■ Ref. [b|] only explored the space of the parameters which 
have a direct effect on the low-multipole temperature spectrum, 2 namely the primordial 
parameters A s , n s , k* and the cosmological constant (which can be re-expressed in our 
basis as Qa = 1 — u m /h 2 ). Other parameters were kept fixed at uj m = 0.135, uj^ = 0.022, 



^ur grids contain the points Lu b SE {0.019, 0.022, 0.025, 0.028}, u cdm € {0.08, 0.10, 0.12, 0.15, 0.18}, 
h £ {0.5, 0.6, 0.7, 0.8, 0.9}, r g {0.0, 0.1, 0.2, 0.3, 0.4}, n s £ {0.80, 0.86, 0.92, 0.98, 1.04, 1.10, 1.16, 1.22, 
1.28} and for the cut-off scale fc* G {0, 10" 3 , 2 x 10" 3 , 10 x 10~ 3 } 

2 The parameters Ha and n s affect both small and large multipoles. So, for consistency, the temperature 
power spectrum should be fitted with no restriction on r, u> m and cot, in order to have any possibility of 
compensating for the effects of (Oa, n 3 ) at large I. 
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r = 0.17. These constraints correspond to the best-fit values (with a running n s ) for the 
"WMAPext+2dFGRS" data set which includes data from other CMB experiments and 
large scale structure measurements. A priori, there is no justification for imposing such 
constraints on a smaller data set (WMAP TT alone). For instance, we find that in the 
absence of a cutoff (k* = 0), the maximum likelihood in the full six-dimensional parameter 
space has an effective x 2 °f 972 for 893 degrees of freedom, while imposing the previous 
constraints raises the best x 2 to 979. When polarization and large scale structure data 
are included in the analysis, fig. 5 shows that the evidence for a cutoff spectrum gets even 
weaker, while values of fc* bigger than 5 x 10~ 4 Mpc -1 are now excluded at 2a. 
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Figure 5: Likelihood of the wave number cutoff fc* for WMAP TT, TT+TE and TT+TE+2dF 
data sets. The TT curve has a slightly higher preferred fc* relative to the other two. The horizontal 
lines refer to the la (68 %) and 2a (95 %) confidence levels. 



To see how parameters are correlated, in fig. 6 we also show confidence regions in the 
r-k* plane for the three different data sets, and the maximum likelihood value of each 
parameter for fixed values k*. Using TT alone (fig. 6a), the preference for nonzero k* is 
only slightly bigger than a la effect. When TE polarization data is included (fig. 6b), the 
significance is further reduced, and the anticipated need for larger optical depth values is 
manifested. By examining the most likely values of the other parameters as a function of 

fig. 6d, we see that an increase in k* and r can be compensated by an increase in n s , 
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h, Mb, and a decrease in u c & m . Inclusion of large scale structure data (fig. 6c) breaks this 
degeneracy, but does not much improve the low statistical significance of the determination 
of k*. 




Figure 6: (a)-(c): 1 and 2 a confidence regions in the t-K plane for WMAP TT, TT+TE and 
TT+TE+2dF data sets, (d) Most probable values of other parameters as a function of fc*, in the 
case of the TT+TE+2dF data sets. 



To give a more precise idea of the statistical significance, we tabulate the best fit 
parameters and x 2 °f the fit for the three cases in Table 1. As expected, the TT data by 
themselves show the strongest preference for a cutoff, but the change in x 2 1S still only 
2.6. When the TE data are included A% 2 falls to 1.8. With the inclusion of 2dF data, 
A^ 2 = 2.0. Thus the model with no cutoff is within lAa of the minimum x 2 point. 



The above results were obtained using the ansatz (3T) for P(k). We repeated the 
analysis for primordial spectra of the forms shown in figs. 2a and 2b, as well as a sharp 
step- function cutoff. They give no better nor worse a fit relative to (3.1). 
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x 2 
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1429.5 


1458.7 


1456.7 


d.o.f. 


1342 


893 


892 


1342 


1341 


1374 


1373 



Table 1: Maximum likelihood parameters for the three different data sets, with and without cutoff 
fc*. WMAP best fit is from table 1 of @. 



4. Summary and conclusions 

We have investigated whether the anomalously small power of large-angle CMB temper- 
ature anisotropics is indicative of new physics that could suppress the primoridal power 
spectrum of inflaton fluctuations at low wave number by introducing a cutoff k*. We ana- 
lyzed the complete WMAP data set including polarization, with and without the addition 
of the 2dF galaxy power spectrum data. Allowing for a variation of the six standard cos- 
mological parameters characterizing a spatially flat universe, we find a marginal preference 
for a nonvanishing cutoff scale of k* = (2.7 ± 1) x 10 -4 Mpc -1 , in agreement with ref. ||12|| . 
Like them we find that the likelihood is nongaussian for low k (fig. 5) so that k* is bounded 
to be above zero at a confidence level of only 59 — 64 % (depending upon whether 2dF 
data is included). By contrast, ref. |]l^] finds A;* > at the 87 — 93 % c.l. Moreover their 
most likely value of /c* is larger than ours, (4.9 ± 1) x 10 -4 Mpc -1 . We have explained the 
reasons for the differences between these results and our own. 

Thus we conclude that at present the motivation from the low quadrupole and octopole 
moments for a power spectrum with an infrared cut-off is quite weak. Nevertheless we have 
pointed out correlations between the parameters which will be relevant for fitting the low 
quadrupole with such models. In particular there is a tendency for optical depth r to 
increase in order that the low-/ polarization data remain consistent with the model in the 
presence of a low-A: cutoff. 

How can we reconcile the small statistical significance of the need for a cutoff with 
the WMAP collaboration's statement that the probability of having such low power at 
large angles is only 2 x 10 -3 ? The answer lies in the highly nonGaussian nature of the 
likelihood function for the C/'s [19, 20]. Expressing the temperature anisotropy as 5T = 
J2im a irnYi m (<j),0), the a/ m 's have a Gaussian distribution, but the multipole moments, 
C( = Sm(l a 'm| 2 ) ' nave the distribution [21] 



dCi ' C, theo ' 1 ' 
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For large I this can be well-approximated by a Gaussian, but for low I it is quite asymmetric 
about the average value x = 1. From this fact we can easily reconcile the small probability 
2 x 10~ 3 with the lack of a significant improvement in the x 2 when considering models with 
cutoff power spectra, as we now explain. 

Let us review the method which was used to arrive at WMAP's small probability, 



as described in section 7 of [15|. From a large ensemble of models in the vicinity of the 
maximum likelihood model, one generates simulations of the data which include the effects 
of cosmic variance and the sky cut which is applied to the actual WMAP observations. 
Among these realizations, one then counts the number whose value of the quadrupole 
moment C2 is less than or equal to the measured C% hs , and compares to the total number 
in the ensemble. If we neglect the experimental noise (which is much smaller than cosmic 
variance at 1=2) and the sky cut (which correlates the quadrupole with other multipoles), 
then for the quadrupole this amounts to computing 

P[C 2 < C 2 obs ] = ^_ Sgfc ca Q.008 (4.2) 



f™dx 



dC 2 



where we used the observed value &2 hs = 123 and the theoretical one C^ 1160 = 1200. Evi- 
dently inclusion of the "experimental complications" decreases this probability somewhat, 
but at least we can roughly understand its small order of magnitude. 

Such an analysis could in principle be repeated including a cutoff power spectrum. We 
can roughly anticipate the results using the same approximation as in equation (4.2), but 
replacing the ACDM integration bound xacdm = 0.1 by ~ 0.2, since the cutoff model 
is able to reduce the theoretical quadrupole by about a factor of 2 (further reduction comes 
at the expense of dragging down the higher multipoles too much). Then we would obtain 
P\p2 < C^ s \ = 0.04. This is a five- fold increase over the low probability in ( |4.2j ), but 
it is still low. Neither model satisfactorily explains the observed low quadrupole, unless 
large statistical fluctuations are invoked. Because of the nongaussianity of the likelihood 
function, only a much more radical suppression of C^ theo could significantly increase this 
probability. 

This explains why in our Bayesian analysis, we do not see a big difference between the 
best x 2 value for fc* 7^ and that for k* = 0. In this approach, we are making a different 
comparision, namely: 

A X 2 ~ -2 hn— (xacdm) - In— (x fc jj -1.6 (4.3) 

This crude estimate does not replace the global seven-parameter analysis that we per- 
formed, but it gives a very good approximation for the improvement in x 2 which we 
find. From this argument we conclude that there is no real discrepancy between WMAP's 
low probability 2 x 10 -3 and the seemingly much higher probability found by us and by 
||, [H, 13 1 . However we also conclude that both the standard ACDM model and the low-fc 



cutoff models are rather poor fits to the observed quadrupole, and to do better one should 
find a way to more effectively suppress the theoretical value of the quadrupole moment. 



- 10 - 



Note added: After the first version of this work appeared, ref. [22| proposed a re- 
alization of the kind of cutoff spectrum we have considered here. Their claim of a 2.5a 
signal (in version 1 of their paper) is subject to the same observations as we have made 
with regard to [13], since they apply exactly the same analysis. We also became aware of 
ref. 



\ 24 1, which considered theoretical implication of the low quadrupole moment when 
it was first measured by COBE. 
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